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A confining, symmetry-preserving, Dyson-Schwinger equation treatment of a vector ® vector con¬ 
tact interaction is used to formulate Faddeev equations for the nucleon and A-baryon in which the 
kernel involves dynamical dressed-quark exchange and whose solutions therefore provide momentum- 
dependent Faddeev amplitudes. These solutions are compared with those obtained in the static 
approximation and with a QCD-kindred formulation of the Faddeev kernel. They are also used to 
compute a range of nucleon properties, amongst them: the proton’s cr-term; the large Bjorken-a; 
values of separate ratios of unpolarised and longitudinally-polarised valence u- and d-quark parton 
distribution functions; and the proton’s tensor charges, which enable one to directly determine the 
effect of dressed-quark electric dipole moments (EDMs) on neutron and proton EDMs. 

PACS numbers: 14.20.Dh; 13.40.Gp; 12.38.Lg; ll.lO.St 


I. INTRODUCTION 


A Poincare-covariant Faddeev equation for baryons 
was introduced in Ref. [l|. In principle, it sums all pos¬ 
sible quantum field theoretical exchanges and interac¬ 
tions that can take place between the three dressed- 
quarks that characterise a baryon’s valence-quark con¬ 
tent; and first, rudimentary computations were described 
in Ref. 0. Numerous analyses of baryon properties us¬ 
ing this framework have subsequently appeared, with the 
level of sophistication and breadth of application increas¬ 
ing steadily, e.g. Refs. [§-l^. The computational effort 
required to solve the Faddeev equation has naturally in¬ 
creased as the kernels have become more complex, so that 
one can consume significant resources in formulating and 
tackling even straightforward problems [H, [H, [T^ . 

On the other hand, numerous practical and realistic 
simplifications present themselves. Chief amongst these 
is suggested by the observation that bound-state studies 
which employ realistic quark-quark interactions (^ . 
predict the appearance of nonpointlike colour-antitriplet 

S ark-bquark (diquark) correlations within baryons |26l - 
. Consequently, the baryon bound-state problem may 
be transformed into solving the linear, homogeneous ma¬ 
trix equation depicted in Fig.[TJ The veracity of this 
approximation was established in Ref. [l^ : and it has 
yielded a wide variety of novel predictions for baryon 
structure and interactions [l3. Ha. [T^. [20l - l^ . Notably, 
empirical evidence supporting the presence of diquarks 
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in the proton is accumulating 

The diquarks within baryons are correlated by dressed- 
gluon exchange and their structure is therefore de¬ 
scribed by a Bethe-Salpeter equation (2fil - l^ . These 
correlations are thus fundamentally different from the 
rudimentary, elementary diquarks introduced roughly 
fifty years ago in order to simplify treatment of the 
three-quark bound-state [sil, [s^ . As highlighted by the 
shaded region in Fig.[I] the two-body correlation pre¬ 
dicted by modern Faddeev equation studies is not frozen: 
all dressed-quarks participate in all diquark clusters . 
The baryon spectrum produced by this Faddeev equation 
should therefore possess significant overlap with that of 
the three-quark constituent model and bear no simple 
relationship with that of the quark-felementary-diquark 
picture. 

Given that solving the Faddeev equation can demand 
extensive numerical analysis, many authors have sought 
simplifications. A drastic approach is to replace the 
Bethe-Salpeter kernel for the quark-quark system by 
a momentum-independent interaction, a motivation for 



FIG. 1: Poincare covariant Faddeev equation, is the Fad¬ 
deev amplitnde for a baryon of total momentnm P — Pq+Pd- 
The shaded rectangle demarcates the kernel of the Faddeev 
equation: single line, dressed-quark propagator; F, diquark 
correlation amplitude; and double line, diquark propagator. 
(See Sec.[TT]for details.) 
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which may be found in the Nambu-Jona-Lasinio model 
[H, [13 • On the other hand, following a body of re¬ 
cent work [illliBI it may also be viewed differently. 
Namely, since gluon propagation is characterised by a 
mass-scale ~ 0.5 GeV, which screens infrared field 
modes [H, idU - lbll l , a vector (g) vector contact-interaction 
can be a realistic representation of the quark-quark scat¬ 
tering kernel in QCD so long as the momentum scales 
being resolved are smaller than mg. Indeed, a symmetry¬ 
preserving Dyson-Schwinger equation (DSE) treatment 
of this interaction yields many results that are practi¬ 
cally indistinguishable from those obtained with the most 
sophisticated kernels that have thus far been employed 

[liii]. 

Recall now that quark exchange within the Faddeev 
equation kernel, highlighted by the shading in Fig.[TJ 
guarantees that no single dressed-quark is treated dif¬ 
ferently from the others. This exchange also provides 
additional binding within the baryon. Here, in the con¬ 
text of a contact vector g) vector quark-quark interac¬ 
tion, which produces diquark correlations amplitudes (F 
in Fig.[T]) that are momentum-independent, dynamical 
quark-exchange is the only feature of the Faddeev ker¬ 
nel which can introduce momentum dependence into the 
baryon Faddeev amplitude, 'k. This complication, too, 
may be avoided if one employs a “static approximation” 
for the exchanged quark, viz. [3 


S{P) 


1 

i"f ■ p + M 


1 

M ’ 


( 1 ) 


where M pz 0.4 GeV is the momentum-independent 
dressed-quark mass obtained in solving the gap equation 
using a contact interaction that provides an efficacious 
phenomenology. Using Fq. (ED. or variants thereof, the 
Faddeev equation collapses to a simple algebraic expres¬ 
sion, whose solution provides both a baryon’s mass and 
a momentum-independent Faddeev amplitude, which ex¬ 
presses the relative strength of various possible diquark 
correlations in the bound-system [s^ - 

Today, Eq. (P) is often a convenience, not a necessity. 
It is employed by some authors because it can assist 
with the development of intuition about the bound three- 
valence-quark system and its properties. Furthermore, 
where comparisons with more sophisticated studies can 
be made, such as for the nucleon, A-baryon, and the 
Roper resonance [T^ lUj, H^l, the static-approximation 
provides an equally good description of baryon masses. 
On the other hand, it provides a flawed description of 
hadron elastic and transition form factors because it leads 
to extreme hardening (^ . and fails completely in de¬ 
scribing the internal structure of the Roper resonance 

M- 

Herein we therefore analyse the nucleon and, to a lesser 
extent, A-baryon Faddeev equations using the contact- 
interaction but eschewing the static-approximation, with 
an aim of determining whether the increased complexity 
brings benefits which outweigh the loss of simplicity. We 
formulate the associated Faddeev equations in Sec.lllland 


describe a procedure for their solution in Sec. nni The 
results for the masses and amplitudes are presented and 
discussed in Sec. lYl 

As noted at the outset, the Faddeev equation is impor¬ 
tant because its solutions enable computation of numer¬ 
ous observable properties of baryons. Hence, in Sec. m 
we present results for the nucleon cr-term, Bjorken-x = 1 
values for the separate ratios of unpolarised and longi¬ 
tudinally polarised u- and d-quark parton distribution 
functions (PDFs), the proton’s tensor charges, and de¬ 
rived results for the neutron and proton electric dipole 
moments (EDMs). We conclude in Sec. lVIl 


II. CONTACT-INTERACTION AND THE 
EADDEEV EQUATION 

A. Gap Equation 

The basic pieces of the Faddeev equation in Fig.[T] 
are the dressed-quark and -diquark propagators, and the 
diquark Bethe-Salpeter amplitudes. The dressed-quark 
propagator is obtained from a gap equation: 

5-i(p)= (* 7 -p + m‘’“) + E(p), (2a) 

S(p) = 9^Dg_^{p - q)^-fgS{q)^T^{q,p). (2b) 

Here is the gluon propagator; F;^, the quark-gluon 
vertex; a symbol representing a Poincare invariant 
regularisation of the four-dimensional integral, with A 
the regularisation mass-scale; and is the current- 

quark bare mass. We employ a confining, symmetry¬ 
preserving DSE treatment of a vector (g) vector contact- 
interaction, which is implemented by writing 

g'^Dg,^{p -q)= , (3) 

mQ 

where ma = 0.5 GeV is a gluon mass-scale typical of 
QGD and the fitted parameter am = 0.36 tt is com¬ 
mensurate with contemporary estimates of the zero- 
momentum value of a running-coupling in QCD [^ . 
This interaction is embedded in a rainbow-ladder (RL) 
truncation of the DSEs [H, [111 , viz. one writes 

^v{p,q) = lu (4) 

in the gap equation and in the subsequent construction of 
the Bethe-Salpeter kernels. Notably, with realistic inter¬ 
actions [ 23 , the RL truncation is known to be reliable 
for properties of the following ground-states: isospin- 
nonzero-pseudoscalar- and vector-mesons, and octet and 
decuplet baryons [stI - I^ . 

Using Eqs. dSD, O, the gap equation becomes 

S~\p) =il-p + m+^^ f jg.S{q)-fg. (5) 

O mQ Jg^g 
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TABLE I: Computed dressed-quark properties, required as 
input for the Bethe-Salpeter and Faddeev equations, and 
derived scalar and axial-vector diquark properties. All re¬ 
sults obtained with am = O.SGtt and (in GeV) Air = 0.24, 
Auv = 0.905. N.B. These parameters are taken from the 
spectrum calculation of Ref. . (Dimensioned quantities are 
listed in GeV.) 


Mo 

mu 

Mu 

''^ppO+ 


EqqO 

FqqO 

Eqql 

0.36 

0.007 

0.37 

0.78 

1.06 

2.74 

0.31 

1.30 


where C'^{oj)/uj = C^'iuj) = r(-l,a;r2^) - r(-l,a;ri^), 
with r(Q;, 2 /) being the incomplete gamma-function. So¬ 
lutions are listed in TableHl N.B. It is a feature of Eq. CHD 
that in the chiral limit, m = mg = 0, a nonzero solution 
for Mg := limmj:->.g Mf is obtained so long as aiR, exceeds 
a minimum value. With Aji-^uv as specified in the Table, 
that value is ~ O.IGtt. This appearance of “mass 
from nothing” expresses the phenomenon of dynamical 
chiral symmetry breaking (DCSB), which is the source 
of more than 98% of the mass of visible material in the 
Universe [t^. 


We have assumed isospin symmetry, so that m„ = m = 
rrid and strong interactions do not distinguish between u 
and d-quarks, and use a Euclidean metric, which is de¬ 
tailed in Appendix A of Ref. [l3] . Equation dS]) possesses 
a quadratic divergence; but if that is regularised in a 
Poincare covariant manner, the solution is 

S{p)~^ = i'j ■ p + M, (6) 


where M is momentum-independent and determined by 


M = m + M 


4aiR 
37rm^ JO 


f 


ds i 


1 


M2 


( 7 ) 


We implement a confining regularisation of the 
contact-interaction by following Ref. and writing 


1 


M2 




dre 


-t(s+M^) _ 


S-kM2 


( 8 ) 
, ( 9 ) 


where are, respectively, infrared and ultraviolet reg¬ 
ulators. This will be our definition of Evidently, a 
finite value of tr =: l/Ajj. in Eq. (jH]) implements confine¬ 
ment by ensuring the absence of quark production thresh¬ 
olds in colour singlet amplitudes. It is worth adding to 
this remark. Contemporary theory predicts that both 
quarks and gluons acquire mass distributions, which are 
large at infrared momenta [H, . The gen¬ 

eration of these mass distributions leads to the emer¬ 
gence of a length-scale c ~ 0.5 fm, whose existence is ev¬ 
ident in all modern studies of dressed-gluon and -quark 
propagators and which signals a marked change in their 
analytic properties. In this realisation, which has been 
canvassed by numerous authors {e.g. Refs. and ci¬ 

tations thereof), confinement is a dynamical process that 
may be realised via Eq. @ . 

The interaction in Eq. (I3|) does not define a renormalis- 
able theory and hence A^v := l/^uv cannot be removed. 
Instead, it plays a dynamical role, setting the scale of all 
dimensioned quantities. Using Eq. the gap equation 
becomes 


B. Diquark Correlations 

It has long been known [2^ that scalar and axial- 
vector diquark correlations are overwhelmingly dominant 
in the nucleon because all are positive-parity states and 
the correlation masses satisfy mggg+, TOqqi+ < toat. No¬ 
tably, the A-baryon involves only axial-vector diquarks 
because it is impossible to build an isospin-3/2 state from 
a dressed-quark and an isospin-zero scalar diquark. 

A contact interaction treatment of diquark correlations 
is detailed in Sec. 2.2 of Ref. M- In this case, the scalar 
and axial-vector diquarks are described by the following 
amplitudes: 

r°9(7’) = il5Eqqo{P) + ■ PFqgoiP), (Ha) 

^Flq^{P) = ^llEqq,{P), (11b) 

where = 0. These amplitudes do not depend on the 
quark-quark relative momentum; and in all applications 
they are canonically normalised (g+ = q + P,Sp {q+) := 
dS{q+)/dP^): 

P^ = 2tr r ro; {-P)Sp^ (g+)ro; iP)S{q) , (12a) 

J dq 

P, = |trrrj;j-P)5p,(g+)rj;jP)5(g) . (12b) 

■J Jdq 

The computed values of the correlation mass-scales and 
the constants which characterise the amplitudes are listed 
in TableHl 


C. Faddeev Equation 

Numerous details relating to the Faddeev equation 
treatment of the nucleon and A-baryon are provided else¬ 
where (e.g. Sec. 3 and Appendix C in Ref. [l^). Here we 
merely recapitulate the main ideas, as they are relevant 
to our analysis. 

The nucleon Faddeev amplitude is 

-H ^f2-k T-a , (13) 


M = m + M 


4aiR 

iiririQ 


C‘"(M2) , 


( 10 ) 


where the subscript identifies the bystander quark and, 
e-5- ^ 1,2 are obtained from d'g by a cyclic permutation 
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of all quark labels. As remarked above, the nucleon is 
composed from a sum of scalar and axial-vector diquark 
correlations: 

4'3 {jpj ,aj,Tj)= Aff + , (14) 

with {pj , aj , Tj) representing the momentum, spin and 
flavour labels of the quarks constituting the bound state, 
and P = Pi + P 2 + P 3 is the system’s total momentum. 
The scalar diquark piece in Eq. (11^1) is 


where Viyp is a 4 x 4 Dirac matrix, which describes the 
relative quark-pseudovector-diquark momentum correla¬ 
tion within the A'*'’’', and Un(P) is a Rarita-Schwinger 
spinor (see Appendix A of Ref. [ij for details). 

The Dirac-matrix structure of Eqs. (Ha, (USD may be 
expressed as follows: 

2 

S{l;P) = J2sn{l;P)T’^{l;P), (24a) 

n—1 


Ar° (K,a„r,)= T^ (K) {K) [S{1; P)u{P)] 


^3 


(15) 


where 5 is a 4 x 4 Dirac matrix, which describes the rela¬ 
tive quark-scalar-diquark momentum correlation within 
the nucleon, the spinor satisfies 


u{P){i'-f ■ P + rriN) = 0 = {i'f ■ P + mN)u{P) , (16) 

with rriN the nucleon mass, K = pi + p 2 ='■ P{i 2 }j ^ ■= 
{-P{i 2 } + 2p3), and 


A°'^(X) 


1 


-I- m 


2 

qqO+ 


(17) 


is the scalar diquark propagator. The colour antisymme¬ 
try of 'I'3 is implicit in T*^ , with the Levi-Civita tensor, 
Ecicaca, expressed via the antisymmetric Gell-Mann ma¬ 
trices, viz. eciC2C3 = when 


8 

M;(;;P) = ^aU;;P)T;(/;P), i = +,0, (24b) 

n—3 

where Sn{l\ P) and a\{l] P) are scalar functions of P, l-P, 
and the 4x4 matrices {P,T^\i = l,2;j = 3,..., 8} 
are defined in Eq. (lAll) of Appendix!^ The A-baryon 
amplitudes may be expressed similarly: 


8 

v,p{i-p) = Y,ui-p)T:p{i-p), ( 25 ) 

n—1 

where fn{l, P) are scalar functions and the 4x4 matrices 
{rip I j = 1,..., 8} are defined in Eq. (IA2I) . 

One can now express the Faddeev equation for 4/3: 


S{k;P)u{P) 

JAp{k;P)u{P) 


{H^ = i\\ = iX^} . 

The axial-vector component in Eq. (1141) is 

■AAg ^ (K, Qfi, Ti) 

4- 1 ^^1^2 4- 

= frl{K) Al,iK)[Al{l;P)u{P)] 


(18) 


r3 

Oi 3 


(19) 


where Ai, is a 4 x 4 Dirac matrix, which describes the rela¬ 
tive quark-pseudovector-diquark momentum correlation 
within the nucleon, the symmetric isospin-triplet matri- 


= -4 


Mik,l-P) 


I dl 


S{l;P)u{P) 

Aiil;P)uiP) 


where 


Mik,l;P) = 


Moo {Moi)i 

{MioYp {MiiYY^ 


with 


(26) 


(27) 


ces are 


0 

0 

= T°\lqq)S^{lqq-kq) 


-1- 

II 

0 

+ 

CO 

li 

II 

0 

1 

CO 

(20) 

{Moi)i 

Xr°\-kqq)S{lq)A°\lqq), 

= t=rl\lqq)S^ilqq-kq) 

(28a) 

and, with T^,,(Ar) = 6 ^ 1 , + KpK^^/rUqq^^ , 



Xr°\-kqq)Silq)AlUlqqY 

(28b) 

, TiiK). 

^ ^ ^qql+ 

(21) 

{MioYp. 

= T°^ilqq)S^{lqq-kq)P 
Xrl\-kqq)Silq)A°-^ilqq), 

(28c) 

The A-baryon contains only an axial-vector diquark: 

(MiiYY. 

= t=rl\lqq)S^ilqq-kq)P 



(22) 


Xrl\-kqq)Silq)Altilqq). 

(28d) 


and when computing the mass and Faddeev amplitude 
one can focus on the A++, owing to isospin symmetry, 
in which case: 


Vl = 




t+T^l 


AlliK)[V,pil-,P)up{P)] 

Oil Q!2 


T3 

0:3 


( 23 ) 


where: Ig = I + P/3, kg = k + P/3, Igg = —I + 2P/3, 
kqq = —k -I- 2P/3; f = r(P)"'’ C, with C = 7274 being 

the charge conjugation matrix, = —C; and the super¬ 
script “T” denotes a transposing of all matrix indices. 

At this point one can use isospin symmetry to define 
A{k-,P) := A°(fc;P) = —M+(fc; P)/-\/2 and therewith 
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simplify Eq. (1261) : 


S{k;P)u{P) 

Afj.{k; P) u{P) 

S{l;P)u{P) 

AAl;P)u{P) 

where, with all entries referring to i, j = 0, 


= - 4 / M{k,l-P) 

Jdl 


M{k,l-P) = 


■Moo ^{Moi)v 


The A-baryon Faddeev equation is 
Vxp{k]P)up{P) 

pA 


= 4 


Mtp{k,l;P)Vp,il-,P)u^{P), 


Idl 


with 


where the form of the entries is obvious, given the de¬ 
scription above and Eq. (IA6I) . e.g. 


(29) 


XZ" = - kg)T°^i-kgg) 

X S{lg)A°^{lgg)Tnil;P)A+{P) ■ (36) 

Following a similar procedure, an analogous equation 
for the A-baryon is readily obtained: 


(30) 


</)^(fc; P) = 4 r ^ P), (37) 

Jdl _ _1 


( 31 ) 


where 


= trfZ{k;P)rl\i,,)s^{i,, - k,)fi^{-k,,) 


p K'^qq)^ Vqq ^qJ^X \ '^qqJ 

X S{QAfMr:Al-.P)Pt{P) ■ (38) 


B. Regularised Equations 


M%{kAP) = t+T];{i,,) 

X - k,)t+fx{-k,,)S{k)^]Jp{k,) . (32) 

III. DEFINING AND SOLVING THE FADDEEV 
EQUATIONS 


As is typical for bound-state equations founded on a 
contact interaction, Eqs. (ESI, EB involve divergences, 
which we tame by using the confining regularisation pro¬ 
cedure described in connection with the gap equation, 
Eqs. El) ®- Foi' illustration, consider the first entry on 
the right-hand-side (rhs) in the equation for (/)(^(fc;P), 
viz. 


A. Scalar Equations 


Since we are not going to employ the so-called static- 
approximation, Eq. the Faddeev equations must be 
solved numerically, with momentum-dependent results 
for all the functions in Eos. (l24l) . (1251) . To begin that 
process, we define a set of Dirac-matrix-valued projec¬ 
tion operators, Eqs. (IA3I) . (IA4I) in Appendix in order 
to convert the linear, homogeneous, Dirac-matrix-valued 
Faddeev equations into a set of equations for these scalar 
functions. Left-multiplying Eq. with the projection 
operators, right-multiplying with ■u(P), and forming the 
spinor trace, one obtains: 

4>Z{k-P) = -A rY.^T<k^{hP), (33) 

Jdl n=l 


where each kernel entry is a function of P, k-P, 
I ■ P, k ■ I, and we have defined 


K{k;P) := 


rt^rnn . 

Aiv •— 


_ J 

^ U 77 ) 

(fc;P), m 
,{k;P), m 

= 1,2 
= 3,.. 

.,8 ’ 

(34) 

n^mn 

A .00 ’ 

m = 1, 2 , 

n = 1, 

2 


OAmn 

A.01 5 

m = 1, 2 , 

co" 

II 

e 

...,8 

(35) 

OAmn 

A.10 5 

m = 3,... 

00 

= 1,2 

OAmn 

, A.11 5 

m = 3,... 

00 

II 

00 



f x]^^nk.p)= f 

Jdl Jd, 

where, using Eq. ([M)) . 


A r\ni 

D 1 D 2 D 3 


!M^\k,l;P)^ 

9 i {k,P) , 


= trTl(fc)r° (lqq)\-i{n -Igg-l - kq) -H My 
X {-kqq){-ij ■ Iq -|-M)ti (/; P) A+(P) , 


D^=ll + M\ 

D2 = {Iqq — kq)"^ M'^ , 

DS. = Iqq + TOggO+ • 


(39) 


(40a) 

(40b) 

(40c) 

(40d) 


The explicit form for is cumbersome, so we do not 
include it here; but it is worth detailing our treatment of 
the denominator product: 


1 

D 1 D 2 D 0 



dxdyx— , 


(41) 


where 


P = (1 - x)Di -I- x{l - y)D 2 + xyDs 

= {I + xyk — {x — r]{l + xy))P)'^ + oj , (42a) 

w = M^[l — a;(l — y)] + xy[k^ — 2k ■ P(1 — y — x) 

- (2r]k ■ P + k'^)xy] + x[ml^o+ (1 - y) 

- - [2 - 7 j]r]y + x[l - . 


(42b) 
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Assuming a Poincare-covariant regularisation proce¬ 
dure, one may shift the integration variable as follows 
- xyk^i -I- (cc - 77(1 -I- xy))P^, so that 

= 9{}^{k, I — xyk -I- (x — 77(1 -|- xy))P; P) (43) 


and hence 

oe^kj-p) 


I dl 


= 2 






/ dxdyx 
lo Jdl 


!H.^\k,l;P) 

(/2+w)3 


(44) 


N.B. After the change-of-variables, I ■ P = 0; and in the 
computation of (fc, 1; P) one may then follow Ref. [s^ 
and set: 


— 5 , /2 

'^qq ' "^qq'i ^qq 


^ggl ^gg/i 


- p I 

? J //I 


—^ -P 

X rs II. • 


3 j ^qq^i ^ 

(45) 

At this point we use Eqs. (I8|), dH) and infer 


1 




i-iy 


£:"+! " nl duj 

in which case Eq. (j44l) becomes 
cA am 


Idl 


9C^{k,l-P) ^2 [ dxdyx 

D 1 D 2 D 3 ^ ^ ^ 

pOO 

/ ^AC\kJ;P)‘4yf+cd)cl>y{k-,P), 

Jdl 


(46) 


(47) 


where = / (i'^l/(27r)^ represents the unbounded four¬ 
dimensional momentum integral. The integral on the rhs 
is now free from ultraviolet divergences and exhibits no 
quark production thresholds, i.e. confinement is realised. 


C. Solving the Regularised Equations 

Eull implementation of the procedures illustrated 
above yields a collection of well-defined linear, homoge¬ 
neous integral equations, one set for the nucleon Eaddeev 
amplitude and another set for the A-baryon. These equa¬ 
tions only have solutions at discrete, separated values of 
p2. As usual, therefore, we consider a modified equation, 
which takes the form, using the nucleon as an example: 

A(p2)0(^(fc; P) = -4 ^ KT € (I-, P) ■ (48) 

'’dl 

Equation (H51) has at least one solution for every value of 
p2; and Eaddeev equation solutions for physical bound- 
states are obtained at those values of P^ for which 
A(P^) = 1. Any sensible numerical procedure may now 
be used to locate bound-state masses and determine the 
associated Eaddeev amplitudes. 


TABLE II: Computed masses of the nucleon and A-baryon, 
obtained with the contact interaction parameters listed in Ta¬ 
ble H The quantities Qn.a are described in connection with 
Eq. (I49II. (Dimensioned quantities are listed in GeV.) 


Herein 


rriN 

ruA 

Ref. [IJ 

rriN 

ruA 


gjv = 1, gA = 1 ffjv = 1.28, ffA = 1.73 
1.30 1.14 

1.65 _L39_ 

gjv = 1, ffA = 1 ffjv = 1.18, ffA = 1.56 
1.27 1.14 

1.60 1.39 


IV. RESULTS FROM THE EADDEEV 
EQUATIONS 


We solved the Eaddeev equations using the contact- 
interaction parameters described in connection with Ta¬ 
ble n The computed masses are listed in Table m 
The first column of numerical entries compares our re¬ 
sults, obtained with a full treatment of the contact- 
interaction, with those produced using the original static- 
approximation, Eq. O- It provides a surprise, viz. any 
additional attraction introduced by our complete treat¬ 
ment of the quark-exchange kernel is compensated by 
the appearance of additional spin-orbit repulsion, ex¬ 
pressed in the momentum-dependent Eaddeev amplitude 
by the presence of the components S 2 (fc; P), 04 , 6 . 7,8 (^5 P) 
in Eqs. (l24l) . The calculated masses are therefore almost 
unaffected by eliminating the static-approximation. 

The numerical values of the calculated masses are 
shifted roughly 25% above their respective empirical val¬ 
ues. It is appropriate that the Eaddeev equation in Fig.[T] 
should yield masses that are larg er than experiment be¬ 
cause, as explained elsewhere the kernel in Fig.[T] 

omits all those resonant contributions which may be as¬ 
sociated with the meson-baryon final-state interactions 
that are resummed in dynamical coupled channels mod¬ 
els [tsI - IsoII in order to transform a bare-baryon into the 
observed state. Our Eaddeev equation should therefore 
be understood as producing the dressed-quark core of 
the bound-state, not the completely-dressed and hence 
observable object. The problem here is that the shift is 
too large. 

Analysis of the effect of meson-baryon final-state in¬ 
teractions indicates that they typically produce a 15% 
reduction in nucleon and A-baryon quark-core masses. 
As noted in Ref. [s^, a Eaddeev equation kernel capa¬ 
ble of producing more realistic quark-core masses can be 
obtained through a modest modification of the quark ex¬ 
change kernel, which is herein implemented thus: 


S{p) 


1 

■ p + M 


9n,a 

i'y ■ p + M 


(49) 


Using this expedient, one obtains the results in the last 
column of Table HIl where the parameters pn.a were cho- 
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0 1 



(GeV^) 


FIG. 2: Upper panel - Zeroth Chebyshev moment of the 
dominant amplitude in the scalar-diquark component of the 
nucleon’s Faddeev amplitude, si (p; P) in the notation of 
Eq. (I24D . plotted for the two values of qn in Table El 
solid (blue) qn = 1.28; red (dashed) qn ~ 1.0. Lower 
panel - similar image for the dominant amplitude in the 
pseudovector-diquark component of the nucleon’s Faddeev 
amplitude, a 3 {p;P). 


sen in order to obtain values for the baryon masses which 
are consistent with estimates of the respective quark-core 
masses. 

It is important to determine whether the modification 
induced by Eq. (|49l) has any effect on baryon internal 
structure. That question can be answered by peering 
into the Faddeev amplitudes. In the static approxima¬ 
tion, the effect is modest, e.g. including > 1 via the 
value in Table HD increases the scalar-diquark component 
of the nucleon by 3%. The outcome is similar when the 
quark exchange is dynamical, as evident in Fig. [2] with 
the normalisation fixed such that the zeroth Chebyshev 
moment of Si (p; P) is unity at = 0, the = 0 strength 
of the leading pseudovector component drops by < 10% 
when gM = 1 1.28 and the p^-dependence of both 

moments is only affected modestly. 

Figure [2] also highlights a salient difference between 
our results and those obtained using the static approx¬ 
imation, Eq. O; namely, whereas the static approxima¬ 
tion produces a momentum-independent result, our am¬ 
plitude, obtained with dynamical quark exchange, ex¬ 
hibits strong momentum dependence: for p^ > 3m^ the 
leading Chebyshev moments of the dominant scalar and 
pseudovector amplitudes depicted in Fig. [2] fall as 1/p^, 





Ipl (GeV) 


FIG. 3: Top panel - Depiction of the complete space- 
time dependence of the dominant piece in the nu¬ 
cleon’s eight-component Poincare-covariant Faddeev ampli¬ 
tude, si(|p|, cosS), computed herein using the contact inter¬ 
action with dynamical quark exchange. Middle panel ~ Same 
function computed using the QGD-kindred Faddeev equation 
kernel described in Refs. [U, [^, which differs from that used 
herein by using momentum-dependent dressed-quark masses 
and diquark Bethe-Salpeter correlation amplitudes. Bottom 
panel - Difference between these two functions. 


up to Inp^-corrections. Such behaviour is typical of two- 
body systems in quantum field theory; and here we have 
a quark-diquark system, with the diquark described by a 
momentum-independent Bethe-Salpeter correlation am¬ 
plitude. 

This outcome motivates a comparison, depicted in 
Fig-El between the amplitude computed herein and 
that obtained using the QCD-kindred Faddeev equa¬ 
tion kernel described in Refs. [2l|, which success¬ 
fully unifies the description of nucleon, A-baryon and 
Roper resonance properties and differs from ours by us- 




































ing momentum-dependent dressed-quark masses and di¬ 
quark Bethe-Salpeter correlation amplitudes. Whilst the 
spacetime dependence of the dominant piece in the nu¬ 
cleon’s eight-component Poincare-covariant Faddeev am¬ 
plitude, si(|p|,cos0), computed using the contact inter¬ 
action with dynamical quark exchange is quantitatively 
different from that obtained with the more realistic ker¬ 
nel, the functions nevertheless show marked qualitative 
similarities. For example, the support for both is concen¬ 
trated in the forward direction, cos 0 > 0, so that align¬ 
ment of p and P is favoured; and the amplitude peaks 
at (IpI ~ Mat/G, cosd = 1), whereat pq ~ P/2 fv p^ 
and hence the natural relative momentum is zero. In 
the antiparallel direction, cos 9 < 0, the support for both 
these functions is concentrated at \p\ = 0, i.e. pq ~ P/3, 
Pd ~ 2P/3. Figure [3] shows that using the static approx¬ 
imation in the contact-interaction Faddeev equation, in 
which case the analogous plot would simply depict the 
curve si(|p|,cos0) = 1, provides a description of nucleon 
structure that is badly flawed in connection with any 
probe sensitive to the nucleon interior; whereas imple¬ 
mentation of dynamical quark exchange in the contact- 
interaction Faddeev equation yields a significant improve¬ 
ment in the description of a baryon’s internal structure. 
Notwithstanding this, it should be borne in mind that 
when computed using a QCD-kindred kernel, the func¬ 
tions in Fig. [3] fall as l/p'^, up to Inp^-corrections, as one 
would expect of a three valence-body system. 

V. PROBING THE NUCLEON 



FIG. 4: Vertex which ensures a conserved electromagnetic 
current for on-shell baryons described by the Faddeev ampli¬ 
tudes, ^ij, described in Sec. im As in Fig.m the single line 
represents S{p), the dressed-quark propagator. Sec. Ill Al and 
the double line, the diquark propagator, Sec. lII Bl F is the 
diquark Bethe-Salpeter amplitude, Sec. lII Bl and the remain¬ 
ing vertices are described in AppendixlBl- the top-left image 
is Diagram 1; the top-right, Diagram 2; and so on, with the 
bottom-right image, Diagram 4. 

The result in Eq. (l5^ is approximately one-half that 
reported in computations using a QCD-derived interac¬ 
tion [sil, 113 and contemporary simulations of lattice- 
regularised QCD [H, dll. As explained elsewhere [^ . 
this mismatch exposes a defect of the contact interaction, 
viz. it produces rigid, momentum-independent diquark 
Bethe-Salpeter amplitudes, an artefact which leads to a 
weaker m-dependence of the diquark (and hence nucleon) 
masses than is obtained with more realistic kernels. Con¬ 
sequently, Eq. (|5^ is an underestimate of a^- 


A. Sigma Term 


B. Canonical Normalisation 


Following Refs. [8l|, dll, it is straightforward to ob¬ 
tain the nucleon’s scalar charge by using the Feynman- 
Hellmann theorem. This method skirts the need to com¬ 
pute the canonical normalisation constant for the nu¬ 
cleon’s Faddeev amplitude because one need only com¬ 
pute the variation of the nucleon’s mass produced by a 
change in current-quark mass. With the Faddeev equa¬ 
tion kernel constructed as described in connection with 
Eqs. (0), dUl), one obtains 

mAr(m) ™-AMeV 112 -b 2.85 m-b 8.54 (50) 

and hence 


The cr-term is an exception. Typically, one must calcu¬ 
late the normalisation before reporting a physical value 
for any observable. The canonical normalisation constant 
for the Faddeev amplitude associated with an isospin 
multiplet is fixed by ensuring that the zero momentum- 
transfer = 0) value of the electric form factor con¬ 
nected with a charged member of the multiplet is equal 
to the electric charge of that state [s^l. In connection 
with the nucleon, whose Poincare-covariant electromag¬ 
netic current may be written [Q = Pf — P/) 

u{Pf)[-i^F^{Q^) + ^a^,Q.F2{Q^)\u{P .), (53) 


drrifq 

m=7MeV 


2.73. 


(51) 


This result is quoted at an hadronic scale appropriate to 
our formulation, viz. (h = 0.39 ± 0.02 GeV (see Ref. [^ . 
Appendix E). On the other hand, the product 


CTAT = m2.73 = 19MeV (52) 

is independent of scale. (Inclusion of meson-baryon loop 
effects is likely to increase the result in Eq. (1531) by 
roughly 7MeV [U.) 


this means imposing Fi{Q^ = 0) = 1 for the proton. 

The nucleon current is detailed in Ref. Q ; and using a 
vector 0 vector contact-interaction, that current may be 
described in terms of the four diagrams depicted in Eig.|4l 
(The so-called seagull terms vanish in our case because 
the diquark Bethe-Salpeter amplitudes are momentum- 
independent.) Diagram 4 vanishes at = 0 and hence 
does not contribute to the normalisation. Its analogue 
does contribute to the tensor charge, however. If one 
chooses to distinguish between quark flavours, the re¬ 
maining diagrams produce fourteen distinct contribu- 
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tions: Diagrams 1 and 2, three each; and Diagram 3, 
eight terms. 

Using the symbols (0,—>-,t) to denote diquark isospin 
labels - 0, scalar [ud] diquark; —>■, pseudovector {ud} 
diquark; and f, pseudovector {uu} diquark, then the dis¬ 
tinct proton diagrams can be expressed as follows: 


^1,2 ^1,2 ^1,2 

^00 ’ 5 

^3 ^3 ^3 ^3 

‘-'00 5 t - O ^ ’ ^Ot ’ ^^0 ’ 

/'*3 /'*3 /'*3 /'♦O 


(54) 


Here, e.g. represents a [ltd] diquark in both the initial 
and final-state proton, with the probe striking a rt-quark; 
and describes a {ud} diquark in the initial-state, 
{wit} in the final state, and the probe striking a it-quark 
exchanged between the two. (The neutron diagrams are 
obtained by exchanging it •<->■ d so that “t” becomes 
and the {itit} diquark is replaced with the {dd} correla¬ 
tion.) Naturally, isospin symmetry reduces the number 
of truly independent computations to just seven: 

^ 1,2 ^ 1,2 ^ 1,2 ^3 
'-'00 ’ ^ ’ ^00 5 

oc oc ; 


and current-conservation makes one of these redundant. 
Our calculation of the normalisation is sketched in Ap- 
oendixIBI 


TABLE III: Selected predictions for the xb = ^ value of the 
indicated quantities. The DSE results are computed using 
the formulae in Eqs. dSSl) - IMll: “realistic” denotes results 
obtained with a sophisticated QCD-kindred Faddeev equa¬ 
tion kernel [IJ]; “contact-S” are contact-interaction results 
obtained using a static approximation, described in connec¬ 
tion with Eq. ©; and “contact-D” are the results obtained 
herein, deriving from a Faddeev equation kernel with dy¬ 
namical dressed-quark exchange. The next four rows are, re¬ 
spectively, results drawn from Refs. [Hil: the row labelled 
-frozen reproduces results from a model in which a non- 
dynamical scalar-diquark is used to describe nucleon struc¬ 
ture, i.e. it is not based on a Faddeev equation and hence 
Diagrams 3 and 4 in Fig. 3] are absent. The last row, labeled 
“pQCD,” expresses predictions made in Refs. [9^, |9g, which 
are based on an SU(6) spin-flavour wave function for the pro¬ 
ton’s valence-quarks and assume helicity conservation in their 
interaction with hard-photons (for reference, 3/7 « 0.43). 



TjVft 

^2 

FP 

^9. 

d 

u 

Ad 

Au 

Au 

u 

Ad 

d 

A/ 

A? 

DSE-realistic [2^ 

0.50 

0.29 

-0.12 

0.67 

-0.29 

0.16 

0.61 

DSE-contact-S 

0.41 

0.18 

-0.07 

0.88 

-0.33 

0.34 

0.88 

DSE-contact-D 

0.38 

0.14 

-0.05 

0.83 

-0.33 

0.43 

0.79 


1 

4 

0 

0 

1 

0 

1 

1 

NJL 

0.43 

0.20 

-0.06 

0.80 

-0.25 

0.35 

0.77 

SU(6) 

2 

3 

1 

2 

1 

4 

2 

3 

1 

3 

0 

5 

9 

CQM 

1 

4 

0 

0 

1 

1 

3 

1 

1 

pQCD 

3 

_ 

1 

5 

1 

5 

1 

1 

1 

1 


C. Valence-quark Distributions 

It is now possible to exploit a connection between the 
= 0 values of elastic form factors, i.e. the Faddeev 
amplitude’s normalisation, and the dimensionless struc¬ 
ture functions of deep inelastic scattering at Biorken-x =: 
xb = 1. As remarked elsewhere [13, EH i whilst all 

familiar parton distribution functions (PDFs) vanish at 
xb = 1, ratios of any two need not; and, under DGLAP 
evolution the value of such a ratio is invariant. 

Thus, e.g. with dv{xB), Uv{xb) the proton’s d, u valence- 
quark PDFs, the value of hma;j 3 _>i dy{xB)/uv{xB) is an 
unambiguous, scale invariant, nonperturbative feature of 
QCD. It is therefore a keen discriminator between frame¬ 
works that claim to explain nucleon structure. Further¬ 
more, Xb = I corresponds strictly to the situation in 
which the invariant mass of the hadronic final state is 
precisely that of the target, viz. elastic scattering. The 
structure functions inferred experimentally on the neigh¬ 
borhood Xb — 1 are therefore determined theoretically 
by the target’s elastic form factors. 

In this connection, consider the current depicted in 
Fig.m Since diquarks are soft, the only contributions 
which survive at xb = 1 are those from Diagrams I and 
3, in which the probe interacts with an isolated quark. 
Each piece from these diagrams appears with a strength 
determined by the proton’s Faddeev amplitude, which 


expresses the effect of weightings derived from both the 
amplitude’s spacetime- and isospin-dependence. These 
properties are expressed in Eq. (jB32|) , which yields the 
following probabilities for finding a quark of flavour / = 

M, d at Xb = 1: 

f^P 

, Pf = /i^u + Ad). (56) 

Inserting computed values for each of the elements in 
Eq. (IB32I) , one finds 

= 0.88, Pd = 0.12, lim^4^=0.14, (57) 

XB—yl Uy[XB) 

which corresponds to P^/Pf = 0.38 at xb = 1. 
These results are collected in Table uni The value of 
dv{xB)/uv{xB) in Eci. (l57)l is 20% smaller than 
that computed using the contact interaction in tandem 
with the static approximation [3^ . In part this is be¬ 
cause the static approximation kills Diagram 3 and that 
diagram favours hard u-quarks; but the reduction oc¬ 
curs mostly because removing the static approximation 
strengthens the Diagram 1 scalar-diquark contribution 
(hard u-quark only) relative to that from the heavier 
pseudovector-diquark. (This can be seen via the discus¬ 
sion of Fig.O) 
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The observations that open this subsection have also 
been exploited in order to deduce a collection of sim¬ 
ple formulae, expressed in terms of the diquark appear¬ 
ance and mixing probabilities, from which one may com¬ 
pute ratios of longitudinal-spin-dependent u- and d-quark 
parton distribution functions on the domain xb ~ 1: 


Pd^) + {Pu^ 

,p _ — Pu^} + {Pd-f — Pd^) 

^ 4(P„^ -I- ) + (Pdf + Pdi) ^ 


(58a) 

(58b) 


where the probabilities for the different quark flavours to 
have helicity aligned (f) or opposite (),) to that of the 
proton are: 

Pu, = P^'^ + (59a) 

Pu, = PSf + (59b) 

Pd, = (59c) 

Pd, = ^PP’^ + \P^'^ ■ (59d) 

The first line of Eq. (|5^ can be understood once one 

understands that is the probability for finding a u- 
quark bystander in association with a scalar [udj-diquark 
correlation in the proton. Owing to Poincare covari¬ 
ance, this term expresses a sum of quark-diquark angu¬ 
lar momentum = 0 and = 1 correlations 

within the nucleon. With = 0, the bystander 

quark carries all the nucleon’s spin. On the other hand, 
the = 1 correlation contributes to both the par¬ 

allel and antiparallel alignment probabilities of the by¬ 
stander quark: 2[u(l\j^u\ud]_^u^ 0 [ud\j^u['ud]_^u^. The rel¬ 
ative strength of these terms is fixed by solving the Fad- 
deev equation and expressed thereafter in the Faddeev 
amplitude: d'o+ ~ 4’l=o + '4’l=i, so that, converting the 
amplitude to probabilities, 


pPfS _ pPiS _|_ pPfS 

Pu\" = ^i=o + 2 iAl=oV'l=i + (60) 

Pi: =Wl=i- 

With the Faddeev equation used herein, PP’'* = 0.82, 
V’L=o = 0.67, V'L=i = 0.24 c/. V'L=o = 0.55, 'ipL=i = 0.22 
in Ref. [2l| and tjjL=o = 0.88, 'ipL=i = 0 in Ref. [s^. 

The other two quantities in Eqs. (15^ are PP>“, pP'”^^ 
which, respectively, gauge the probability that the pho¬ 
ton interacts with an axial-vector diquark component of 
the nucleon or induces a transition between diquark com¬ 
ponents of the incoming and outgoing nucleon. Our dy¬ 
namical dressed-quark exchange Faddeev kernel gener¬ 
ates PP’“ = 0.18 and PP’™ Ri 0; and we list results for 
numerous ratios in Table IIIIl (As remarked above, the 
result PP’'" R! 0 owes to the generally small magnitude 
of each Diagram 3 term and interference between their 
contributions, see Table [BT]) 


It is worth reiterating that the results in Table |TTT] high¬ 
light that no single ratio is capable of completely distin¬ 
guishing between distinct pictures of nucleon structure. 
Conversely, they show that a comparison between exper¬ 
iment and different predictions for the combination of all 
tabulated quantities provides a very effective means of 
discriminating between competing descriptions da. 


D. Tensor Charges 

With the normalisation computed, one may also read¬ 
ily calculate the proton’s tensor charges, which are de¬ 
fined via {q = u,d): 

{P{p,<^)\qcrp,,q\P{p,a)) = 6Tqu{p,<T)(Tfj,^u{p,a), (61) 

where |P(p, cr)) is a state vector describing a proton with 
momentum p and spin al The derived isoscalar and 
isovector tensor charges are: 

^ = Stu + drd , ^ = Stu — dxd. (62) 

Importantly, the tensor charges are scale-dependent 
quantities, as explained, for instance, in Appendix F of 
Ref. [13. The values decrease uniformly as the resolving 
scale is increased. We compute the results at (h and 
use one-loop evolution equations in order to also report 
values at (2 ■= 2 GeV. 

The nucleon’s tensor interaction is qualitatively iden¬ 
tical to the photon-nucleon interaction depicted in Fig. 2] 
If one distinguishes between quark flavours in this case, 
there are sixteen distinct contributions to each tensor 
charge because Diagram 4 is nonzero. Using isospin sym¬ 
metry, this tally is reduced to only eight. The analysis is 
illustrated via consideration of Diagram 4 in AppendixICl 
Notably, 


= CoZ (63) 

and CZ = 0 because a scalar correlation cannot possess 
a tensor charge. 

Our computed results for the proton’s tensor charge 
at the hadronic scale Ch are summarised in the upper 
panel of Table IlYl It is natural to compare the listed 
values with those in Ref. HO] , obtained using the contact 
interaction and the static approximation, and listed in 
Table C3 therein. To facilitate that comparison, we note 
that Diagram I herein equates to the sum of Diagrams I 
and 2 therein. Diagram 2 herein is the sum of Diagrams 
3 and 4 therein. Diagram 3 is absent in Ref. [13 owing to 
the static approximation, and our Diagram 4 equates to 
the sum of Diagrams 5 and 6 therein. The net results of 


^ In the isospin symmetric limit: := 5tu = S^d, Sl^d : = 

(Syd = where the superscripts denote the hadron in which 

the indicated valence-quark resides. 
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TABLE IV: Upper panel - Proton tensor charges evaluated at 
the model scale: (^h = 0.39 ±0.02 GeV, partitioned according 
to contributions from the diagrams in Fig.[4]and summed to 
provide the complete result. Lower panel - Same quantities 
evaluated using the static approximation. The listed values 
serve as an update of the results in Ref. p3| . 


Dynamic 

Su 

Sd 

(0) 

9t 


Diagram 1 

0.72 

-0.039 

0.69 

0.76 

Diagram 2 

0.14 

0.027 

0.16 

0.11 

Diagram 3 

0.12 

-0.050 

0.075 

0.17 

Diagram 4 

-0.19 

-0.19 

-0.38 

0 

Total 

0.79 

-0.25 

0.54 

1.05 


Static 

Su 

Sd 

(0) 

9t 


Diagram 1 

0.56 

-0.036 

0.53 

0.60 

Diagram 2 

0.29 

0.059 

0.35 

0.23 

Diagram 3 

0 

0 

0 

0 

Diagram 4 

-0.25 

-0.25 

-0.50 

0 

Total 

0.61 

-0.23 

0.38 

0.83 


both calculations are semi-quantitatively similar. How¬ 
ever, there are some line-item discrepancies. We have a 
nonzero result for Diagram 3, which adds to both 5u, 6d, 
whereas this contribution is naturally absent when the 
static approximation is used. There is also a difference 
between the Diagram 2 contributions. This owes to a 
combination of the inadvertent omission in Ref. of a 
factor of two in the probe-diquark vertices, which we have 
restored, with the corrected results listed in the lower 
panel of Table IIVI combined with our use of Eqs. (IB21I) 
in simplifying computation of the currents. 

It is worth analysing the microscopic origin of the pro¬ 
ton’s tensor charges in our calculation. The dominant 
contribution to Spu arises from Diagram 1: tensor probe 
interacting with a dressed u-quark when a scalar diquark 
is the bystander. The next largest piece is produced by 
Diagram 4, in which the tensor probe excites a transi¬ 
tion between the scalar and pseudovector diquarks; but 
this is largely cancelled by the sum of Diagram 2 (tensor 
probe interacting with the pseudovector diquark with a 
dressed-quark spectator) and Diagram 3 (two-loop dia¬ 
grams in which the tensor probe interacts with a dressed- 
quark “in-flight”). Qualitatively equivalent interference 
was seen in Ref. [2^, so that one may still conclude that 
6 tu directly probes the strength of DCSB and hence the 
strong interaction at infrared momenta. Owing to anal¬ 
ogous interference between Diagrams 1-3, Diagram 4 in 
Fig. H] is responsible for the bulk of Srd. 

Notably herein, in contrast to Ref. [2^, Spd ^ 0 even 
in the absence of pseudovector diquark correlations: ow¬ 
ing to Diagram 3, the tensor probe can interact with a 
dressed d-quark exchanged during the breakup and ref¬ 
ormation of a scalar diquark. However, as already noted, 
the contribution is small, viz. this term produces just 


20% of the Diagram 3-dd entry in the upper panel of Ta¬ 
ble IIVI and all other entries would be zero in the absence 
of pseudovector diquarks. The magnitude of bxd may 
therefore still be interpreted as a measure of the strength 
of pseudovector diquark correlations in the proton. 

In ascribing an error to our final result for the ten¬ 
sor charges, we follow the reasoning in Ref. [2^. Namely, 
since one generally finds that systematic treatments of 
the contact interaction yield results for low-momentum- 
transfer observables which are practically indistinguish¬ 
able from those produced by RL studies that employ 
more sophisticated interactions [13, IsH - lil] : and analy¬ 
ses of hadron physics observables using the RL trunca¬ 
tion and one-loop QCD renormalisation-group-improved 
(RGI) kernels for the gap and bound-state equations pro¬ 
duce results that are normally within 15% of the experi¬ 
mental value [1^, we therefore attach a relative error of 
15% to the results in Table IIVI Hence, our predictions 
are: 




M 


5tu 

0.79(12) 


A A 

07 ’a grp gp 

-0.25(4) 0.54(8) 1.05(16) 


(64) 


Given that our computed value of the proton’s cr-term 
is too small by a factor of roughly two, one might be 
concerned by the size of our error assignment. This 
concern can be allayed by first noting that the small 
size of (7p can be tracked directly to an underestimate 
of the diquarks’ cr-terms [l^l- These quantities mea¬ 
sure the rate-of-change of a “charge” (the diquark mass, 
in this case) associated with variations in an external 
source. Although those rates-of-change are underesti¬ 
mated, the masses are not. Thus, whilst our framework 
might produce rates-of-change for the tensor charges that 
are too small, the values of the tensor charges themselves 
should be accurate within the usual error associated with 
rainbow-ladder truncation. 

The results in Eq. (IMl) are quoted at the model scale. 
In order to make a sensible comparison with estimates ob¬ 
tained in modern simulations of lattice-regularised QCD, 
those results must be evolved to ^2 = 2 GeV. We therefore 
list here the results obtained under leading-order evolu¬ 
tion to ((2 = 2 GeV, obtained via multiplication by the 
factor 0.794, explained and computed in Appendix F of 
Ref. [13: 


5tu 


C2 


0.63(9) 


Srd gp^ 
-0.20(3) 0.43(6) 


9t^ 

0.83(12) 


(65) 


We use the one-loop expression owing to the simplicity of 
our framework. Employing next-to-leading-order evolu¬ 
tion leads simply to a 25% increase in <Ch with no material 
phenomenological differences. 

The predictions in Eq. (1651) are compared in Fig.[5]with 
phenomenological analyses 1^ . 1^ that benchmark a 
proposed JLab experiment l97| and res ults obtained us¬ 
ing numerous other methods [ToO-lIol. Evidently, our 
predictions match, within er rors, the m ost recent results 
obtained using lattice-QCD [lOOl Il0l| . which appear at 
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FIG. 5: Flavour separation of the proton’s tensor charge: “1” 
- illustration of anticipated accuracy in planned JLab ex¬ 
periment [^, with central values based on the analysis in 
Ref. [ 9 ^; “2” - results drawn from Ref. [ 9 ^; “3” phenomeno¬ 
logical estimate in Ref. [9^ “4” - prediction herein, Eq. (I65II : 
“5” - corrected results from Ref. [ 2 ^, drawn from the lower 
panel of Table nYi and evolved to (" 2 ; “6-14” - estimates from 
Refs. [l0Cll4l08l| . respectively. The bands drawn from “4”-“7” 
are described in connection with Eq. IMl)- By way of con¬ 
text, we note that were the proton a weakly- inter acting col¬ 
lection of three massive valence-quarks, then [l08l| the quark 
axial and tensor charges are identical, so that 5tu = 4/3 and 
= —1/3 at the model scale. These values are located at 
“15”. If one assumes that C,h may reasonably be assigned as 
the scale of that model, then those values evolve to Stu = 1.06 
and Srd = —0.26, which we also plot at “15”. 


“ 6 ” and “7” in Fig.jSJ A weighted combination of the 
most recent DSE- and lattice-QCD results yields the fol¬ 
lowing estimates, drawn in Fig. [5] 

= 0.70 ± 0.03 , (5Td=-0.21 ±0.01. ( 66 ) 

Another interesting point is highlighted by a compar¬ 
ison between our predictions and the values obtained 
when the proton is considered to be a weakly-interacting 
collection of three massive valence-quarks describ ed by 
an SU(4)-symmetric spin-flavour wave function [l08|: 

= 2e„ and = Cd cf. our results, Eq. (IMl) . 

Stu = 0.59(2e„), Sxd = 0.75(ed). It is thus apparent 
that the presence of diquark correlations in the proton 
Faddeev amplitude materially suppresses the magnitude 
of the tensor charge associated with each valence quark 
whilst simultaneously increasing the ratio STd/Sxu by 
approximately 30%. 


E. Electric Dipole Moments 

In typical exte nsions of the Standard Model, quarks 
acquire an EDM |l09L Ill0l| , i.e. an interaction with the 


photon that proceeds via a current of the form: 

975 0 -^ 1 /<?, (67) 


where dq is the quark’s EDM and here we consider 
q = u,d. The EDM of a proton containing quarks which 
interact in this way is defined as follows: 

{P{P, cr)\J^^^\PiP, cr)) = dp uip, a)-f5ap^u{p, a), ( 68 ) 
where 


= duu(x)j5af,^u{x) + ddd{x)-f5ap^d{x). (69) 
At this point, using a simple Dirac-matrix identity: 

75*^^!^ — 2^/ii^cK/5^Q/3 ; (70) 


one can write 


yyEDM 


2 ^ 


; ± dd da add 


(71) 


It follows that 


{P{p,a)\J^,^^\P{p,a)) 


du 5 tu + dd Srd 


uip, cr)75 uip, a); (72) 


namely [ 2 ^, the quark-EDM contribution to a proton’s 
EDM is completely determined once the proton’s tensor 
charges are known: 

dp = du Stu + dd Sxd . (73) 


With emerging experimental techniques, it is possible to 
place competitive upper-limits on the proton’s EDM us¬ 
ing storage rings in whi ch polarized particles are exposed 
to an electric field 111 |. 

An analogous result for the neutron is readily inferred. 
In the limit of isospin symmetry. 


(Nip, a)\uap,,u\Nip, a)) = (P(p, cr)|dcr^i.d|P(p, cr)), 

{Nip,a)\dap,,d\Nip,a)) = {Pip,a)\uaf,uu\Pip,a)); 

(74) 


and hence dn = du St d + dd Stu. Using the results in 
Eq. ((Ml) , we therefore have 

dn = -0.25 du + 0.79 dd, dp = 0.79 - 0.25 dd ■ (75) 

It is worth contrasting Eqs. dza with the results one 
would obtain by assuming that the nucleon is simply a 
collection of three massive valence-quarks described by 
an SU(4)-symmetric spin-flavour wave function. Then, 
by analogy with magnetic moment computations, a pro¬ 
cedure also made valid by Eq. (|70)) ; 

dn — 3 du ± 3 dd , dp — 2 du 3 dd , (76) 

values which are roughly 50% larger than ours. 
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The impact of our tensor-charge predictions on 
beyond-Standard-Model phenomenology may be eluci - 
dated, e.g. by following the analysis in Refs. [mini. 
In this connection it is worth remarking that the pos¬ 
sibility of a s-quark contribution produces cons iderable 
uncertainty in estimates of nucle on E DMs |ll4| . largely 
because its size is very uncertain 115l |. Therefore, even a 
rudimentary DSE estimate of this contribution could be 
useful. Such may be obtained via a simplified treatment 
of meson-loop corr ections to the quark gap equations, 
as used elsewhere [HE [m to estimate the proton’s 
strangeness-magnetic-moment and -cr-term. Following 
that reasoning, one is led to (5 ts(C2) ~ 0.02^^ = 0.009. 


VI. CONCLUSION 

We employed a confining, symmetry-preserving, 
Dyson-Schwinger equation treatment of a vector ® vector 
contact interaction in order to formulate and solve Fad- 
deev equations for the nucleon and A-baryon in which 
the kernel involves dynamical dressed-quark exchange. 
These are the first contact-interaction calculations to pro¬ 
duce momentum-dependent Faddeev amplitudes. Pre¬ 
vious contact-interaction studies have imposed a sup¬ 
plementary condition on the Faddeev equation kernel, 
viz. a “static approximation” [Fq. ([T|)], which leads to 
momentum-independent amplitudes. 

So far as computed masses of the nucleon and A- 
baryon are concerned, eliminating the static approxi¬ 
mation has little effect [Table HI] . On the other hand, 
the impact on the Faddeev amplitudes is dramatic. In 
stark contrast to static approximation results, contact- 
interaction dynamical quark exchange produces ampli¬ 
tudes that compare far more favourably with those ob¬ 
tained using kernels built from elements that possess 
momentum-dependence typical of QCD [Fig. [3]. This 
marked improvement can potentially lead to a better de¬ 
scription of those baryon properties which are sensitive 
to momentum scales that exceed the dressed-quark mass, 
M ^ 0.4 GeV. That can be tested, e.g., in a study of the 
nucleon’s first radial excitation and the related transition 
form factors, for which the contact-interaction-|-static- 
approximation yields poor results. 

Regarding the proton’s cr-term, the contact interac¬ 
tion result is one-half that obtained using more realistic 
interactions [Fq. (l5^ ]. whether or not the static approx¬ 
imation is employed. The small value is thus a defect of 
the contact interaction itself, which may be traced to the 
rigid nature of the diquark masses and Bethe-Salpeter 
amplitudes obtained with this interaction. 

Implementing dynamical quark exchange in the 
contact-interaction Faddeev equation kernels yields re¬ 
sults for the large Bjorken-x values of the separate ra¬ 
tios of unpolarised and longitudinally-polarised valence 
u- and d-quark parton distribution functions that are no¬ 
ticeably different from those obtained using the static ap¬ 
proximation. The differences owe primarily to an increase 


in strength for the scalar diquark component of the nu¬ 
cleon. Our dynamical results [Eqs. (l57l) ] should be viewed 
as a truer representation of the contact-interaction’s pre¬ 
dictions. 

The increase in strength for the scalar-diquark com¬ 
ponent of the nucleon, induced by dynamical quark ex¬ 
change, also affects the computed values of the proton’s 
tensor charges, dyd [reported in Fq. (1351) and com¬ 
pared with other estimates in Fig. [5]: it acts primarily to 
increase 5tu. Our calculations confirm that the presence 
of diquark correlations in the proton reduces the size of 
6tu, StiI compared with results obtained in simple quark 
models whilst simultaneously increasing \STd/dTu\; ver¬ 
ify that 6tu is a direct measure of the strength of DCSB 
in the Standard Model; and emphasise that \STd/STu\ di¬ 
minishes with pP’°-/pP’‘^^ i.e. the ratio of pseudovector- 
and scalar-diquark interaction probabilities. 

With this analysis we have completed the first im¬ 
provement promised in Ref. [IE and so must address the 
question posed in the Introduction, viz. Does the in¬ 
creased complexity which accompanies dynamical quark 
exchange in the Faddeev kernel outweigh the loss of sim¬ 
plicity inherent in using the static approximation? If 
one is driving toward a realistic picture of a wide range 
of hadron physics observables, including those sensitive 
to probe momenta greater than the dressed-quark mass, 
then the answer must be affirmative. This question has 
a natural extension, however: Is there merit in contin¬ 
uing to use a contact interaction when computational 
resources are beginning to enable the use of realistic in¬ 
teractions in the study of baryons? Here the answer is: 
Yes, depending upon the problem in hand; but the char¬ 
acter of those problems is rapidly evolving. 

Finally, returning to the proton’s tensor charges, the 
next step should be computation using the approaches of 
Refs. ™ order to obtain continuum predictions 

with a direct connection to QCD. 
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Appendix A: Collected Formulae 

The matrices in Eos. (l24)l . (|25|l . which express the 
Dirac-matrix structure of the positive-energy nucleon and 
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A, are 

T^(;;P) = I4x4 , t'^{ 1;P) = i'y ■ , 

: Tt.Q‘\P) = ■ ^"^5 , 

{l;P) = -iP^ilb , tI{1-P) = P^7 • [-^75 , (Al) 

T'®(^;-P) = ■ It - , 

where + 1 ■ P = Iv + 1 ■ P Pv, P = 1, 

= —1; and 

Tlp{l]P) =Sp^, 

Tlp{h = ^ ( 27^-4 - ^ 

tIpP, P) = -iVSPpl^j ■ It , 

T^p{l-,P) = V3Ppi^, (A2) 

Ttp{l]P) 

TupP,P) = , 

rlpP, P) = ^ + 3/^/^ , 

We using the following projection operators in arriving 
at the scalar-valued integral equations in Sec. lIIIl - for the 
nucleon, 

fi(/;P) = il4x4, fpl-,P) = -^j-l^, 

^pP^P) = , t4(/;P) = -^ 757^7 . f-J-, 

Tp{l;P) = -^nbPf., f®(;;P) =-i75P^7-f-^, (A3) 

^®(^; P) = -^ 75 ( 7 ^ 7 • It - ip); 

and for the A-baryon, 

rlp{l-,P) 

P.pih P) = -^ (27.^^p^ - 3^-p7 • l^) , 

T^pi^P) = -i^Ppl-t-fl-^, 

f^pillP) = -^Ppi^. (A4) 

Tlp{i]P) = 

Tlp{l\P) = \lplv ^ 

^IpP,P) = -\iui ■ ■ 

Tlp{i\P) = + ipi^ - ^ipi^i■ 


These projectors are defined such that, for the nucleon: 

tr [P{1)t\1)A+{P)] = i,j = 1,2, (A5a) 

tr[r;(0r^(0A+(P)] *,j = 3,...,8, (A5b) 

tr[P(0T^(0A+(P)] =0, z = l,2,j = 3,...,8, (A5c) 

where the positive-energy projector is defined via 

2A/A+(P) := u{P,s)u{P,s) = (-iy • P + M). 

S —± 

(A6) 

In connection with the A-baryon, 

Tr [PpAlpf^pmUP)] = = 1,. •. ,8, (A7) 

where 

PUp) = 

- hpA. + iPpPu - k^P^.1. - P.7m)]A+(P). (A8) 

Appendix B: Nucleon-Photon Vertex 

1. Diagram 1: Normalisation 

When a symmetry-preserving regularisation of the con¬ 
tact interaction is employed, one finds 

Tl{Q^ = 0)=lp, (Bl) 

i.e. the dressed-photon-quark vertex preserves its bare 
form at zero momentum-transfer . Consequently, the 
Cqq (scalar-diquark) contribution to the nucleon normal¬ 
isation constant, 9^, can be written 

— 

QloA+iPhpA+{P)= %,A+iP) Sik-,-P) 

Jdk 

X 5(fc,)e„7M5(fc,)A0+(A:,,)5(fc;P)A+(P), (B2) 

where S{k\P) is the [ud] scalar-diquark component 
of the nucleon’s Faddeev amplitude in Eq. (I24al) and 
S{k;—P) = C'^S{—k;—P)C] S{kq) is the dressed-quark 
propagator described in Sec. Ill Al and = (2/3). 

After some algebra, following the pattern in Sec. 1111 Bl 
Eq. (IB2I) yields 

pi poc 

Qjo = 2e„fA4 / dxx ‘E 2 {k'^ + ^o)‘T>ooih P), (B3) 

Jo J dk 

where, with n{tj = —zy • (. + M, 

T’mik;P) = tri7^A+(P) 

X S{k- -P)n{kg)-fpn{kg)S{k; P)A+(P), (B4a) 
k = k + {^ — x)P , (B4b) 

Wo = xM'^ + (1 - x)rrP^^Q+ - a:(l - x)m1p . (B4c) 
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Inserting our gjv = 1-28 computed nucleon mass and 
Faddeev amplitude into Eq. (IB3I) , we obtain 

Q^oo = 0-0163 eu% =: Q°eu% - (B5) 


Diagram 1 also represents two other cases, viz. a 
dressed-quark struck with either {ud} or {uu} axial- 
vector diquarks as the spectator, Cf^ respectively. 

Consider first the case of a w-quark struck with a {ud} 
spectator. Analysis of the type described above yields 


Qn = 



where 



(E2(fc2+wi)®i\(fc;P), 


(B6) 


= tri7^A+(P)A„(fc;-P) 

X n{kq)-f^n{kq)T^pikqq)A/ 3 {k;P)A+{P), (B7a) 

wi = xM"^ + (1 — ^)'^'^qqi+ ~ 2^(1 ~ x)m'jq . (B7b) 

Inserting our gjq = 1.28 nucleon mass and Faddeev 
amplitude into Eq. (IB6I) . we obtain 

= 0.00148 eu% =: . (B8) 

Owing to isospin symmetry, the other case, d-quark 
struck with a {mm} spectator, contributes (e^ = 

-e„/2) 


{ed/eu)Qlii-V2f = -Qli. (B9) 

Plainly, therefore, diagrams with axial-vector diquark 
spectators do not contribute to the proton’s normali¬ 
sation. More generally, in fact, they contribute noth¬ 
ing to the proton’s electromagnetic form factors [s^: 
+ Cl-f = 0. 


2. Diagram 2: Normalisation 

In this instance, the first contribution we consider is 
a M-quark spectator to a photon-scalar-diquark interac¬ 
tion, 

— 

QloA+{Phf.A+{P)= %A+{P) S{k;-P) 

Jdk 

X A°\kqq)V^{kqq)A°\kqq)S{kq)S{k-P)A+{P), (BIO) 

where, owing to a Ward-Green-Takahashi identity main¬ 
tained by a symmetry-preserving regularisation of the 
contact interaction [37|, the dressed photon-scalar- 
diquark vertex at zero momentum transfer is {Cud = 

Cu -f Cfz) 


^d^kqq) — 2 Cud kqq ^ . (Bll) 

Using this result and the now standard procedures, 
Eq. (|B10I) yields 


(B12) 


Diagram 2 also represents two other cases, viz. a 
dressed-quark spectator to a photon interacting with ei¬ 
ther a {ud} or {uu} axial-vector diquark, re¬ 

spectively. The u-quark spectator contribution is; 


/•A 

Ql,A+{P)j^A+{P)= %A+{P) / Au,ik-,-P) All, {kqq) 


Jdk 

X V^,c'fi>{kqq)Al%{kqq)S{kq)A{k-,P)A+{P) , (B13) 


which involves the dressed photon-pseudovector-diquark 
vertex at zero momentum transfer [37| : 

U^,a/3(A;gg) = ‘2 Cudkqq ^T^^{kqq) . (B14) 

Using this result, Eq. (IB 131) yields 

Ql, = 0.000879 Cud %=:'D^Cud %■ (B15) 

Owing to isospin symmetry, the analogous result for 
the d-quark spectator diagram, is (e„„ = 2e„ = 4e„d) 


(e„„/e„d)(-V^)^gn = 82)^e„d % . (B16) 


3. Diagram 3: Normalisation 


If quark flavours are distinguished, then this image in 
Fig. m corresponds to eight two-loop diagrams. The sim¬ 
plest example involves a scalar-diquark in the initial and 
final states, so that a dressed d-quark is struck “in-flight” 
by the photon, Coq- 

QIoA+{P)-/^A+{P) = %A+{P)edgj,[ [ S{t,-P) 

Jdk Jdl 

X S{kq)A°^ {kqq)T°'" {lqq)[S{-kx)ldS{-kx)r 
X f°\-kqq)S{lq)A°^ {lqq)Sil-,P)A+{P) , (B17) 

where kx = k + I — P/3. After some algebra, Eq. (IB 171) 
can be recast in the following form: 


rA rA 


Ooo “ ^dgx^ 


A? 


-00 


dk Jdl -D1D2D3D4D5 


(B18) 


with 


Dl — k/^p M^, D2 — /Cqq -I- , 

D3 = 1^ + M^, D4 = , (B19) 

D5 = k^ P , 

id 

= tlil^A+iP)Sik-,-P)S{kq)A°\kqq)T°\lqq) 

X [ni-kx)-fM-kx)rf°^{-kqq) 

X n{lq)A°\lqq)S{l-,P)A+{P). 


Qla = 0.0152 Cud % =■ 'D^Cud % ■ 


(B20) 
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In the algebraic evaluation of , the following replace¬ 
ments are used sequentially: 


k 

k 


2 ^ 777 ^ 

qq ^ "‘ggO+’ 

-P 

99^^ 3 


/2 _^2 
tgg "‘‘qqO+1 

1 Ip 

'■99 A* ^ 3 A*' 


(B21) 


As usual, the next step is to combine the denomi¬ 
nators using a Feynman parametrisation, in which case 
Eq. (IB18I) becomes 

/»A /»A 

Qoo = / 

Jdk Jdl 

X 24: f dxidx2dx3dx/^x\x\x\x4^ ^°r°n3 ’ (B22) 

where the variable transformations k ^ k, I ^ I are 
constructed so that Hgg = f{k'^,P,P'^ = —m\). At 
this point one can implement the confining regularisation 
prescription and thereby obtain 


/»oo /»oo 

Qoo ~ / 24 / dxidx2dx3dxix\x\x\x4^ 

J dk J dl Jo 


X ^QoikJ;P) J . 


(B23) 


The precise forms for the elements in Eq. (IB23I) are 
lengthy so we do not present them here. Notwithstanding 
that, their computation is straightforward; and inserting 
our g]\[ = 1.28 nucleon mass and Faddeev amplitudes 
into the expressions one derives, the following numerical 
result is obtained: 


Qlo = 5.34 X 10-^edg%% = 8.75 x 10-^ed% (B24) 

=: X°°ed% . (B25) 


There are seven more contributions; but, owing 
to isospin symmetry, as described in connection with 
Eq. (|55L only two additional computations are necessary: 
{ud} pseudovector-diquark breakup, d-quark struck in¬ 
flight, [ud] scalar diquark recombination {[ud]d:y{ud})^ 
Co_>; and {ud} pseudovector-diquark breakup, d-quark 
struck in-flight, {ud} pseudovector diquark recombina¬ 
tion ({ud} dg,{ud}), The results are: 

= 1.79 X lO-^edlA^ =: X°^ed% , (B26a) 


[ud] Uj{uu} = (7q| = —V2euX^^9{p , (B26b) 

{ud} dg,[ud] = tAo = edX^^9{p , (B26c) 

{uu} Ug,[ud] = (7|o = —V2euX^^ 9{p ; (B26d) 

and 

= -8.8 X IQ-^edlA^ =: X^^ed% , (B27a) 

{ud} Ug,{uu} = = —\/2euX^^ , (B27b) 

{uu} Ug,{ud} = = — V2euX^^ 9{p . (B27c) 


TABLE B.l: Breakdown by diagram of contributions to the 
nucleon’s canonical normalisation. Column 2: Results scaled 
as described in Sec. lB 51 N.B. In relation to Eqs. (IB28II . (IB32I) : 
Diagram 1 corresponds to summing Q,terms; Diagram 2, sum¬ 
mation of ® terms; and Diagram 3, summation of X terms. 
(Actual value of each entry is obtained via division by 10®.) 



i/a:^ 

1 /A^ 

Diagram 1 

10.9 

10.5 

Diagram 2 

7.70 

8.46 

Diagram 3 

-0.0729 

0.157 

Diagram 4 

0 

0 

Total 

18.5 

19.2 


4. Collected Results: Normalisation 

For ease of reference, we gather all independent results 
computed in Secs. IB llIB 31 into Eg. (IB28p . wherein each 
entry should be divided by 10®: 

qO q1 p, 0 rpl j(-00 _;^00 

16.3 1.48 15.2 0.879 0.0875 0.0179 -8.80 x ’ 

(B28) 

The associated contributions to the nucleon’s canonical 
normalisation constant are listed in Table iBTl 


5. Current Conservation 

In a symmetry preserving treatment of the contact in¬ 
teraction, using integration by parts and changes of vari¬ 
ables, one can establish the following Ward identities be¬ 
tween the contributions described above: 

Cl° = ‘D° + 2X°° -b 2(1 - V2)X°^ , (B29a) 

Q} = + 2X°^ + 2{1 - V2)X^^ , (B29b) 

In the static approximation, Eqs. O, (HU), the exchange 
diagrams vanish and Eqs. (IB29I) reproduce the identities 
used in Refs. [H, [H, l44j |. 

At the same level, the neutron should be neutral. Con¬ 
sider, therefore, that in the isospin symmetric limit, the 
neutron’s charge can be computed from the diagrams de¬ 
tailed in Secs. IB IMB 31 bv making the exchange u O d: 

Cn = %[Q?ed + Q}{ed + 2e„) 

+ Pk' S-ud + ‘P^iSud + 2ecid) 

+ X^^Cu + 2X^~^(eu — y/2ed) 

+ A^^(e,-2V2ed)] =0, (B30) 

if Eqs. (IbMI) are satisfied. Numerically, however, they 
are not. That is because the last step in our regularisa¬ 
tion procedure involves the introduction of hard infrared 
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and ultraviolet cutoffs, which naturally spoils the condi¬ 
tions necessary to prove Eqs. (IB29|) . Following Ref. [^ . 
a remedy is straightforward. Namely, we use Eq. (IB29cll 
as a definition of and introduce rescaling factors Tq i, 
with ^0 = 0,0(1-Hro), 2)0 ^ 2)0 = 2)0(l-ro) and 

O^ ^ O^ = 0.^(1 + n), ^ = 2)^1 - n), such 

that Eqs. (IB29I) are satisfied when expressed in terms of 
^0.1, ^0,1^ procedure yields 

ro = -0.0298, n = -0.203, (B31) 


Analysing and combining these integrals in the now 
customary way, then at the hadronic scale C,h = 0.39 ± 
0.02 GeV one obtains 

= -3.68 X = -0.192 . (C4) 

Both C^q contribute half their strength equally to 
Stu, Srd, so that 


and entails ^ = 119X0^ =: rxX°^. N.B. 

These same scaling factors are applied to the analogous 
diagrams that appear in computing the proton’s tensor 
charges. 

The nucleon Faddeev amplitude can now be nor¬ 
malised, following the procedure described in Sec. lVBI 
One has 

Cp = + Serf) 

+ + 2)^ (find + 2e 

+ X°°ed + 2X°^{ed-V2eu) 

-bX^^(ed - 2y2e„)] = 0.0192 (B32) 

^!A^ = 52.2. (B33) 

The contribution from each diagram-number is listed sep¬ 
arately in Table IB. II 


Appendix C: Nucleon-Tensor Vertex 


When a symmetry-preserving regularisation of the con¬ 
tact interaction is employed, one finds [2^ 






(Cl) 


i.e. the tensor vertex is not dressed, a result which owes 
to the inability of any symmetry-preserving regularisa¬ 
tion of a contact interaction to support nonzero relative 
momentum in a quark-antiquark system. 

Following our analysis of the proton’s electromagnetic 
current at = 0, only one type of diagram remains to be 
explicated, viz. Diagram 4 in Fig.|4l In connection with a 
zero-momentum tensor probe, this image translates into 
two terms: with the latter corresponding to 


%\A+{P)a^,A+{P) = %A+{P) / Aa{k-,-P) 

Jdk 

^ ^a/3ikqql+)'^l3fJ,L>ikqql+ , kqq0+)^ ikqq0+) 

x5(fc,)5(fc;P)A+(P), (C2) 

where the tensor transition vertex is {Nf = 2, = 2): 


rA 

I3fiu(kqq\+ ,kqq0+) ~ ~‘2AffN^ tr / S{q p kqqi+) 

J dq 


’ fiiy 


X S{q + kqq 0 +)T°^{kqq 0 +)S{q)rl i-kqqi+). (C3) 


c4 

Ot 


u = S^d = %% 


(C5) 


The complete list of contributions from Diagrams 1-3 
is: 


= 0.0147 (1 + ro) = 0.744, (C6a) 

= -4.86 X 10"^ %{l + ri) = -0.0195, (C6b) 

= 2PX^ , (C6c) 


o 

II 

O 

(C7a) 

= 8.48 X 10-'‘ %{l- ri) = 0.0544, 

(C7b) 

(7"2 _ 

(C7c) 


rr3 _ 
-^00 — 

-1.93 X 10-'‘ % = -0.00992 , 

(C8a) 

( 7-3 _ 

8.8 X 10”^ = 0.00459, 

(C8b) 

II 


(C8c) 

rr3 _ 
-^^0 — 


(C8d) 

II 

ts 


(C8e) 

p^^ = 

-7.87 X 10”® %rx = -0.0487, 

(C8f) 

T4t = 

-V2Pl^ , 

(C8g) 

Tf4 = 

-V2P^^ . 

(C8h) 


In terms of the contributions in Eqs. (IC5l) - (IC8p . the pro¬ 
ton’s tensor charges are: 


Stu = [Pol + 24 ^] 

+ [K4^ + 'z^t] 


+ [^t ■*" 

+ 4-0 + ‘44] + I 

(C9a) 

Srd = P^\ + ^ 24 ^ 



+ [Pol + %4 

+ 24o -b 24 ^] -b 24 o . 

(C9b) 


These expressions yield the results in the upper panel of 
Table IVl 
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